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Abstract. 

We investigate the dynamical behavior of entanglement in a system made by 
two solid-state emitters, as two quantum dots, embedded in two separated micro- 
cavities. In these solid-state systems, in addition to the coupling with the cavity 
mode, the emitter is coupled to a continuum of leaky modes providing additional losses 
and it is also subject to a phonon- induced pure dephasing mechanism. We model 
this physical configuration as a multipartite system composed by two independent 
parts each containing a qubit embedded in a single-mode cavity, exposed to cavity 
losses, spontaneous emission and pure dephasing. We study the time evolution 
of entanglement of this multipartite open system finally applying this theoretical 
framework to the case of currently available solid-state quantum dots in micro-cavities. 
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1. Introduction 

Cavity quantum electrodynamics (CQED) deals with the interaction among photons 
confined in a reflective micro-cavity and atoms or other particles. When a two- 
level system (for example, a two-level atom) is strongly coupled to a cavity mode 
PQ, it is possible to realize important quantum information processing tasks, such as 
controlled coherent coupling and entanglement of distinguishable quantum systems 
[2 El H]. In this respect solid-state devices, and in particular semiconductor 
quantum dots (QDs), utilized as "artificial atoms", are one of the most promising 
architectures for the possibility of miniaturization, electrical injection, control and 
scalability. Recently, thanks to impressive progress in the technology of solid-state 
microcavities, substantial advances have been made towards these goals. The strong 
coupling regime has been reached for the excitonic transition of quantum dots (QDs) 
El [71 [H [9] , and nanocrystals [10] coupled to optical semiconductor cavities, as well 
as for superconducting qubits coupled to microwave cavities [TT]. In all of these 
systems, the cavity-mode quality factor (Q) can be very large while solid-state emitters 
are intrinsically coupled to the matrix they are embedded in. In fact, decoherence 
and phase relaxation unavoidably broaden any transition between the discrete states 
of such artificial atoms. High experimental performances are required to realize 
quantum processors and it is thus important to establish how long a sufficient degree of 
entanglement can be maintained in spite of losses, decoherence and noise. Implications 
are the possibility to store entangled states in solid-state memories and entanglement 
preservation during local operations in quantum algorithms [12]. Understanding how 
the entanglement is transferred from, e.g., a pair of independent initially entangled 
qubits to reservoirs has motivated several contributions in recent years [131 E]- The 
aspect that has mostly drawn attention is the possibility of a complete disappearance 
of the entanglement between the qubits at finite times [15]. The occurrence of this 
phenomenon, termed entanglement sudden death (ESD) [16] and of entanglement 
revivals [17] have been shown in quantum optics experiments [HI [19]. Entanglement 
transfer from atoms to cavity modes leading to entanglement revivals has also been 
studied [20]. The effects on the dynamics of entanglement of cavity losses, spontaneous 
emission and pure dephasing for two-qubit systems have been extensively investigated 
[2T1 [22]. The general problem of the dynamics of entanglement in the simultaneous 
presence of more than one noise has been also studied finding that, for composite 
systems, the additivity of decay rates of single noises is not maintained [23]. The 
main aim of this paper is to investigate the dynamical behavior of entanglement for a 
realistic implementation for quantum computation of two qubits in separated cavities 
where several sources of noise are present. 

Here we consider the entanglement dynamics in a system composed by two initially 
correlated solid-state emitters each strongly coupled to a lossy cavity interacting with its 
reservoir. We include the unavoidable losses due to spontaneous emission into external 
electromagnetic modes distinct from the lossy single-mode of the cavity. These, always 
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Figure 1. (Color online) Schematic representation of the four-partite system (qA, ca\ 
Qb, c B ) 

. The two qubits qA and qs are initially entangled. 

present, additional losses arise from the coupling of the emitter to a continuum of leaky 
modes [23]. We also include pure dephasing noise which plays a relevant role in solid 
state quantum emitters. For example, a QD interacts with the phonons of the matrix it is 
embedded in, giving rise to sidebands in addition to the so-called zero-phonon line (ZPL). 
At sufficiently low temperature yet, the emission in the ZPL remains predominant, 
allowing to model these systems as effective two-level systems subject to additional pure 
dephasing [24]. We also study the entanglement transfer from the two-emitter system 
to the cavity modes. The input-output relations for optical cavities [25] show that the 
entanglement between cavity modes can in principle be measured experimentally by 
collecting photons escaping the cavities. Such two-cavity entanglement dynamics could 
thus be exploited to monitor the entanglement dynamics of cavity-embedded solid state 
emitters. 

The paper is structured as follows. In Sec. |2] we present and solve a model for the 
physical configuration described above. In Sec. [3] we explore the dynamical behavior of 
entanglement for quite general values of the physical parameters of the system. In Sec. H] 
we specialize to the case of currently available solid-state quantum dots in micro- cavities 
giving the characteristic lifetimes of entanglement. In Sec. |5]we summarize our results. 

2. Model 

Our system is composed by two noninteracting subsystems (S = A,B), each consisting 
in a qubit (two-level emitter) q$ coupled to a single-mode cavity c$ in turn interacting 
with an external reservoir r$ (see Fig. 1). The Hamiltonian of the total system is thus 
given by the sum of the Hamiltonians of the two noninteracting subsystems 

H t ot = Ha + H B . (1) 

In each subsystem, we distinguish the bipartite system made by the qubit and the 
cavity from the reservoir r. The Hamiltonian H$ of each part S = A, B reads like (we 
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omit index S) 

H = H Q + H r + H = ^u a z + u c a) a + g(a) <r_ + a <r+) , (2) 

where if Q describes, in the rotating wave approximation (RWA), the qubit-cavity system, 
g is the coupling constant between qubit and cavity, er z denotes the usual diagonal Pauli 
matrix, a± are the two-level raising and lowering operators, a and a) the annihilation 
and creation operators for the cavity mode. We observe that this Hamiltonian model is 
valid for values of g smaller than ojq ~ u c [26]. H r describes the external environment 
responsible of the different noise sources which affect the qubit-cavity system and if; 
the interaction of the latter with the environment. The realistic conditions present in a 
system composed by quantum dots embedded in microcavities are modeled as three noise 
sources: cavity losses, qubit spontaneous emission and pure dephasing mechanisms. 

In the usual master equation approach, considering the Markov approximation and 
an infinite number of bath oscillators, we can describe the dynamics of the qubit-cavity 
system by 

j t P = i[p,H } + £p. (3) 

The Markovian processes are described by the Liouvillian term, that consist of 3 parts: 

C = C cav + C SE + C D ■ (4) 

C cav describes the cavity losses of photon in the reservoir modes and is expressed by the 
following form [12] : 



£>cav = -y(2apa* — a)ap — pa) a) , (5) 

being 7 C the photon escape rate from the cavity to free space. In addition, the qubit 
is subjected to decay via spontaneous emission and losses of coherence. Cse describes 
spontaneous emission in the leaky modes: all the available light modes except the cavity 
one [12], 

T 

£se = y(2a_po- + - a + a_p - pa + a_) . (6) 

Cca V and Cse are based on Hamiltonians in RWA and thus are not valid for arbitrarily 
large values of the decay rates 7 C and 7 q [26]. Finally, 

C D = ^j{a z pa z - p) (7) 



describes pure dephasing processes [T2]. Although the three noises are treated in the 
Markovian limit, the reduced dynamics of the qubit-cavity system can present non- 
Markovian behavior depending on the strength of the coupling constant g with respect 
to the various decay rates. We shall comment quantitatively on this point after Eq. ( TTTj) . 

The optical cavity is an open quantum system, cavity photons can escape it and 
propagate into free space on along an optical fibre until they eventually reach another 
distant quantum system or can be detected. The quasimode approach is able to describe 
in a direct way such a behavior. A relationship between the external fields and the 
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intracavity field may be obtained f25\ 127] in the limit of continuous spectrum. If the 
coupling constant k(u) between cavity and external bosonic modes is independent of 
frequency over a band of frequencies about the characteristic frequency u c , k(u) ~ k: 



where the operators a out and a in are related to the reservoir operators and 7 C = 2itk, 2 . 
Throughout this paper we will consider the case of no input photons, hence, once known 
the quantum state for the cavity mode, it will be possible to calculate expectation values 
and correlation functions for output photons that can be measured experimentally or 
be used as input to transmit the entanglement to distant quantum systems. 

2.1. Procedure 

Being the two subsystems noninteracting, they evolve independently so that we can 
analyze the dynamics of only one subsystem and in turn obtaining the evolution of 
the total four-partite open system [T7J. We will consider initial states with zero cavity 
photons and at most one excitation in each qubit. Eq. ([3]) allows us to compute the 
joint evolution of the total four-partite open system starting from an arbitrary initial 
state where only one excitation is present in each atom. From the knowledge of the 
evolved density matrix, it will be possible to investigate the various reduced dynamics 
of the total system, for example that of the two qubits or of the two cavities. In the 
following we start showing how to compute for each part the time dependent density 
matrix elements, which in general may be different for the two subsystems. 

2.1.1. Dynamics of subsystems Here we consider the dynamics of a single subsystem 
S. We choose the standard product basis B = {|1) = |l q )|l c ), |2) = |l q )|0 c ), |3) = 
|0 q )|l c ), |4) = |0 q )|0 c )}, where |0 9 ) (|0 C )) and (|l c )) are the lower and upper state 
of the qubit (cavity). The dynamics of qubit qs under the effect of the master equation 
of Eq. (j3J) is described by the reduced density matrix 



The time dependent coefficients Pt and pt can be obtained analytically, however the 
presence of pure dephasing gives rise to very cumbersome and lengthy equations. In 
this section we present analytical results only for the case 7d = 0. Numerical results in 
presence of pure dephasing shall be included in next section. Analogously, for the cavity 
modes, the dynamics of the reduced density matrix can be expressed as 



where a zero-photon initial state (pii 6 (0) = Pk^O) = and p 2 2 6 (0) = 1) nas been 
considered. The dynamics of the reduced density matrices is obtained using the 
master equation expressed by Eq. Q, that is appropriate when the reservoir is at zero 




(8) 




(9) 




(10) 
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temperature, the coupling between the cavity and the external modes of the reservoir 
has a flat spectrum in the range of involved frequencies and the qubit is resonant with 
the cavity [281 EH]- We will limit our investigation to this physical condition. In the 
absence of pure dephasing, from Eq. fl3]) one obtains, 



Pt 



Qt 



(Tc+Tq) + 
4 1 



(7c+7q) t 
4 r 



cos(f2t) + 



7c -7q 



sin (fit) 



— sin(fW) 



with P t 

Q = 



p 2 and Q t = q% and where we have introduced the characteristic frequency 



g 2 — ((7c — 7 q )/4) 2 . We point out that the function Q t is not directly connected 
to the decay of the cavity excited state but is linked to the exchange of the initial 
qubit excitation between the qubit itself, the cavity and the reservoir. In this sense, the 
fact that Q t goes to zero when g — reflects its dependence on the initial conditions 
and on the internal qubit-cavity dynamics. Non-Markovian features in the dynamics of 
the subsystems occur for values of g large enough to make f2 real, but such as not to 
compromise the validity of the RWA in our model. 

We notice that these direct relationships between the functions appearing in the 
diagonal and non-diagonal elements of Eqs. fl9]) and f fTOj) fail for 7^ ^ 0. Following 
Ref . [17] , the knowledge of two any single parts dynamics permits to obtain the dynamics 
of the corresponding bipartite system. 



3. Entanglement evolution 

After obtaining in the previous section all the relevant dynamical coefficients, in this 
section we give the explicit expressions of concurrence for some couples of parties 
of the four-partite system. In particular we consider the two-qubit and two-cavity 
entanglement separately for two different initial configurations. The qubits are initially 
in one and two-excitation Bell-like states, while cavities are in their vacuum state. 

We shall restrict our analysis of entanglement dynamics to the two-excitation 
entangled initial states 

|¥)= (a|00) 9 + /3|ll) 9 )|00) c = (a|4) 9 + /3|l) 9 )|4) c , (12) 

where in each ket the first entry denotes a (q or c) state of subsystem A, while the 
second entry a state of subsystem B. Generalization of the results for other (eventually 
mixed) initial states is straightforward. 

From the evolved state \^>t) one finds the reduced density matrices of the bipartite 
system of interest tracing over the degrees of freedom of the noninvolved parties. 
We represent the density matrices in the standard computational basis B = {|1) = 
1 11), |2) = |10), |3) = |01), 1 4) = |00)}. In this way, the two-qubit state at time t is, e.g., 
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given by p* Aqg (t) 

( WP t A P t B 
o 



\ a(3*pfpf 



|/3| 2 P/(1 - Pf) 








;i - p t A )p t 





B 



afipfpf 



a 2 + |/3| 2 (l-P^)(l 



Pf) ) 



.(13) 



In the following we will consider the case of two identical subsystems, hence pf = pf . 
An analogous result can be obtained for the two-cavity system density matrix pf ACB {t) 
just replacing p t and P t with q t and Q t respectively. 

The concurrence corresponding to the density matrix of Eq. [13] in the case of 
identical subsystems is found to be [Hj 



a 



max < 0, 2 



-9(c) 
P14 



q(c) q(c) 
P22 P33 



(14) 



We present numerical calculations for a specific two-excitation entangled initial 
state 1^) with a = 0.8. Figure 2 displays the concurrence dynamics of the qubits C% 
and of the cavity modes C£ in absence (2a) and presence (2b) of pure dephasing. Figure 
2 also displays Q t = (at a), providing information on the detectable output photon flux 
(a' out a ont ) = 7c Qt- in both cases 7 C = 0.3g, the other decay rates being fixed as 7 q = 0.3g 
and 7d = in Fig. 2a while 7 q = 0, 7d = 0.3g in Fig. 2b. For these values, non-Markovian 
features appear in each qubit-cavity subsystem dynamics, leading to similar effects in 
the dynamics of the plotted quantities. The figure shows that pure dephasing affects 
heavily the entanglement dynamics increasing the entanglement decay and enabling its 
sudden death. Although phase noise acts directly only on the emitter, owing to the 
strong coupling between the emitter and the cavity, it affects dramatically also the two- 
cavity concurrence. Figure 3a shows the two-emitter concurrence as function of time 
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Figure 2. Q t = (ala) (green solid line), concurrences of the qubits C|, (black dot- 
dashed line) and of the cavity modes C^, (red dashed line) for the initial state |\E') 
with a = 0.8 and with j c — 0.3g as a function of the dimensionlcss quantity gt in 
non-Markovian regime, (a) 7 q = 0.3<?, 7d = 0. (b) 7 q = 0, 7d — 0.3g. 



and of the amount of phase noise Figure 3b displays the two-cavity concurrence C%. 
We used a = l/\/2, 7 C = 0.17 g, 7 q = 0. The detrimental effect of phase noise on both 
the dynamics of the entanglement of emitters and cavities is evident. 
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Figure 3. Concurrences of the qubits C|, (panel a) and of the two spatially separate 
cavity- modes C% (panel b) as a function of the dimensionlcss quantities gt and jd/g 
for the initial state V E') with a = l/v2 in non-Markovian regime for 7 C = 0.17g, 7 q = 



4. Application to quantum dots under realistic conditions 

Here we specialize to the case where the considered system is implemented by using 
currently available quantum dots as quantum emitters, embedded in separated micro- 
cavities [9]. Solid state microcavities with three-dimensional photon confinement, high 
Q, and small volume mode can be realized by fabricating a photonic crystal slab 
structure with a nanocavity composed of one or more missing air holes. The slab 
incorporates a central layer of low density self-assembled InAs quantum dots [5j [9]. 
Another geometry of particular interest is that of micropillar cavities [HI [30] . In these 
systems the fundamental cavity mode can be coupled to and from the outside with 
a very high coupling efficiency. Moreover, they offer interesting perspectives for the 
implementation of quantum information protocols using charged quantum dots [31] . 
Figure 4a displays the concurrence dynamics of the quantum emitters C\ and of 
the cavity modes C| obtained for two independent cavity-embedded quantum dots 
(a = 0.8). We consider typical system parameters for the state-of-art microstructures 
[9]: 7 C = 100 yiteV, 7d = 30/ieV,7 q = 10 /xeV, g = 110 /ieV. The figure shows that the 
two-dot entanglement after a rebirth survives up to about 40 ps. Quantum operations 
based on all-optical implementations can be performed by means of ultrafast pulses. At 
optical frequencies, pulses of 20-100 fs time-width are currently available. Very recently 
structures displaying higher Q values but with a quite low coupling g have been realized 
|30j . Figure 4b displays the concurrence dynamics of the quantum emitters C|, and of 
the cavity modes C£ obtained by using a = 0.8 and parameters describing this novel 
micropillar structure [30]: 7 C = 20 /xeV, 7d = 12/ieV,7 q = 4 /ieV, g = 16 /xeV. In 
this case both the two-dot and two-cavity modes entanglement increase their lifetime. 
Figure 4 also displays the behavior of Q t = (a^a) which is proportional to the detectable 
output photon flux (a out a ou t) = 7c Qt- 

We notice that in the systems considered in this section, dots and cavities 
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frequencies are ~ 1.3 eV 3> g so that the validity limits of RWA are well satisfied. 
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Figure 4. Q t 



(green solid line), concurrences of the quantum dots C|, (black 



dot-dashed line) and of the cavity modes C£ (red dashed line) for the initial state 
with a = 0.8 as a function of the dimcnsionlcss quantity gt in non-Markovian regime. 
Values of other parameters: (panel a) j c = 100 ficV. 7d = 30/zeV,7 q = 10 /xeV, 
g = 110 ficV and (panel b) j c = 20 /ieV, 7d = 12^eV,7 q = 4 ^eV, g = 16 /u,eV. 



5. Conclusions 

We have studied the dynamics of entanglement in a system made by two solid-state 
emitters, as two quantum dots, embedded in two separated micro-cavities. In addition 
to the coupling with cavity mode, the emitter is subject to spontaneous emission, 
due to the coupling with a continuum of leaky modes, and to phonon-induced pure 
dephasing mechanisms. We have modeled this physical system as a multipartite system 
composed by two independent parts each containing a qubit exposed to cavity losses, 
spontaneous emission and pure dephasing. The numerical results presented here for 
arbitrary values of the physical parameters put forward the impact of pure dephasing 
on the entanglement dynamics of the quantum emitters and of the two-cavity modes. 
Experimental information about the latter can be gathered by detection of the collected 
cavity output field. We have finally applied this theoretical framework to the specific 
case of currently available solid-state quantum dots in micro-cavities. 
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